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CN ■ Abstract 

. Hadronization, the process by which energetic quarks evolve into 

hadrons, has been studied phenomenologically for decades. However, 
' httle experimental insight has been gained into the space-time features 

. of this fundamentally non-perturbative process. New experiments at 

>^ [ Jefferson Lab, in combination with HERMES data, will provide signif- 

icant new insights into the phenomena connected with hadron forma- 
^ ■ tion in deep inelastic scattering, such as quark energy loss in-medium, 

^ , gluon emission, and color field restoration. 



1 Introduction 

Hadronization has been a ubiquitous feature of high energy physics experi- 
ments for more than two decades. The subculture that has evolved around 
the study of this process has succeeded in characterizing the topological fea- 
tures of high energy jets, such as the fragmentation functions that describe 
the probabilities of quarks of a given flavor evolving into particular hadrons. 
The string picture is an example of a successful phenomenology that cap- 
tures the important features over a wide range of energies. Few aspects of 
hadronization, however, are presently understood at the most fundamental 
level, since it is a non-perturbative process. 



*This work is supported in part by the U.S. Department of Energy, including DOE 
Contract No. DE-AC05-84ER40150. 
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One important type of experimental data that is generally not easily 
accessible for hadronization is its space-time development. While some phe- 
nomenological models can be used to predict this development, the pro- 
cess takes place on a microscopic distance scale that is not directly ob- 
servable. However, this type of information can be inferred by implant- 
ing the hadronization process into a nucleus, which acts as a spatial filter. 
Fully formed hadrons in the nuclear medium interact via ordinary hadronic 
cross sections; by contrast, propagating quarks do not appear to interact 
in a way that modifies hadron yields significantly. Therefore, much can be 
learned from studying hadron production in deep inelastic scattering kine- 
matics (DIS) from nuclei of varying sizes. The degree of modification of 
fragmentation functions, for instance, can be studied as a function of nuclear 
size A, energy transfer z/, hadron energy fraction z = Eh/v-, and the momen- 
tum component perpendicular to g, which is labelled pr- These studies can 
yield important insights into the fundamental processes driving hadroniza- 
tion, such as gluon emission and the temporal evolution of the truncated 
color field following hard interactions. 

Several previous measurements of semi- inclusive leptoproduction of hadrons 
from nuclei have been performed. The earliest generation of these were car- 
ried out with electron beams at SLACjTj and muon beams at CERN0 and 
FNAL|3]. Prominent features displayed by these efforts are an attenuation 
of hadrons which increases with A and vanishes with increasing v. Follow- 
ing these early efforts, recent data|ljjn] from the HERMES experiment at 
DESYfB] have inspired a new wave of theoretical interest. Unlike the his- 
torical efforts, these new data offer excellent hadron identification and high 
precision, allowing the flavor dependence of hadron formation to be studied. 

Finally, existing and future measurements at Thomas Jefferson National 
Accelerator Facility (Jefferson Lab) will offer unique capabilities to study 
hadron production in DIS. Measurements at the CEBAF Large Acceptance 
Spectrometer (CLAS) at Jefferson Lab[7| will provide data on the widest 
possible range of nuclear target masses at high luminosity. The first data at 5 
GeV were taken in 2003 with sustained luminosities approaching 10^^cm~^s~^, 
while 11 GeV data are planned for the future, accompanied by an order-of- 
magnitude luminosity increase and improved particle identification using the 
upgraded CLAS ("CLAS++"). 

The scientific topics connected with these studies are also important 
themes in other physics communities. Quark propagation through strongly 
interacting media is a foundational element in predicting signatures of the 
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quark gluon plasma at RHIC through the phenomenon of jet quenching jH]; 
in this connection, it has been explored in ultrarelativistic d-A collisionsjn]. 
It has also been investigated in the Drell-Yan process at FNAL in studies of 
quark energy loss in the nuclear medium [T^. 

2 Experimental Tools 

A primary experimental tool for nuclear hadron formation is the hadronic 
multiplicity ratio: 



where N/^^^ and N^^^ denote the number of hadrons and electrons measured 
in DIS kinematics, is the four-momentum transfer, and cj) is the angle be- 
tween the lepton scattering plane and the photon-hadron plane. The primary 
connection to theory has traditionally been through R\j interpreted in a par- 
tonic framework, for particular flavor combinations. Additional tools include 
energy and angular correlations between outgoing hadrons involved in the 
reaction. There are at least two classes that may be of interest: correlations 
between two energetic hadrons formed directly from the hadronization of the 
struck quark, and correlations between the energetic hadrons and very low 
energy protons. The pt distribution is expected to become broader for larger 
nuclei; this transverse momentum broadening is expected to be large enough 
to be directly measurable |lUj. In addition to these quantities, polarization 
observahles such as single spin asymmetries may also play a role. 

A variety of physical properties can be derived from the above observables. 
Hadron formation lengths^, the characteristic distances over which hadrons 
form, can be extracted from R^j. Quark energy loss in passing through the 
nuclear medium may also be accessible, particularly via transverse momen- 
tum broadening, which is also closely connected with a quark-gluon correla- 
tion function. Correlations between final-state particles are expected to shed 
light on the reaction mechanisms, and the nature of the interaction with the 
nuclear medium, whether partonic or hadronic. 

^In some models, referred to as "production length," the distance to produce a "pre- 
hadron" in a two-step picture of hadron formation. 
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A variety of experimental capabilities are required in order to fully ex- 
ploit the observables discussed above. Access to the detailed multivariate 
dependences of R^j requires a large statistical sample; this can be most 
economically accessed at high luminosity. Isolating the energy and angular 
correlations requires large acceptance, since the kinematic variables of the 
particles of interest are not strongly correlated. Good particle identification 
for several hadron species is required to analyze the flavor dependence of 
observables. A large energy range is needed to understand the u dependence 
as well as test the validity of parton model assumptions. Polarization is 
needed to study such quantities as single spin asymmetries. Accessing all 
of these experimental ingredients will require combining information from 
several different facilities. 

3 Model Descriptions and Contact with Ob- 
servables 

Lattice computations can address nonperturbative processes in QCD. How- 
ever, there are several difficulties in applying this technique to hadronization. 
A natural choice to describe the time development of the process is Minkowski 
space, however, this introduces some significant technical obstacles. In ad- 
dition, the relevant space-time volume needed requires a very large lattice 
compared to calculating static properties of hadrons. These calculations may 
become feasible in the 5-10 year time frame. 

Therefore, to make progress on this topic, modeling is needed. Most 
relevant are models possessing close contact to QCD, and which can predict 
a wide variety of observables. A synopsis of model approaches is given in the 
following section. 

3.1 Models 

Historically, a number of models have been developed to describe hadron pro- 
duction in nuclear DIS. While it is not practical to review the older models 
here, a representative approach by Bialas and Chmaj^^ is instructive. In 
this model, the struck quark evolves into a hadron with a characteristic time 
r. The probability that the propagating object is a quark at a given time t 
after it is struck is Pg = e"*/"^ and the corresponding probability that it has 
become a hadron is simply = 1 — Pg. Phenomeno logical cross sections 
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for the quark and hadron interacting with the nuclear medium are postu- 
lated: aq and ah. The probability of the quark or hadron interacting with 
the nuclear medium is taken as proportional to the density of the medium 
for a spherical nucleus, using a standard density parameterization. R\j can 
then be calculated by numerical integration. This model approach was qual- 
itatively successful in describing the v dependence of R\j for the HERMES 
nitrogen datajlj, but failed to describe the z dependence even qualitatively. 

In a more sophisticated approach, the gluon bremsstrahlung modeljT2] 
divides the hadronization process into several stages. As the struck quark 
propagates, an initial stage of uninhibited gluon emission occurs, followed by 
a second stage of reduced gluon emission. The suppression of gluon emission 
in the second stage is imposed to conserve energy, since the propagating quark 
begins with only a finite amount of energy, v. At some point in the second 
stage, a color dipole is formed from the struck quark and the last emitted 
gluon. The color dipole subsequently evolves into a hadron. The role of 
the nuclear medium is to modulate the rate at which the color dipole turns 
into a hadron, including the effect of color transparency; in later versions of 
the model, medium-stimulated gluon emission is included, which has a small 
effect on the prediction for R\j. The assumptions of the model are valid for 
pions with z > 0.5. Its predictions agree well with the HERMES nitrogen 
data in both the z dependence and the u dependence. 

A more recent model approach to describing hadron production in nu- 
clear DIS is based on a continuous series of developments since the late 
1980's addressing jet quenching in relativistic heavy ion reactions|H]. In this 
approach, the prediction for the behavior of R'lj is based solely on medium- 
stimulated gluon emission; hadronization is assumed to occur outside the 
nuclear medium. The calculation of R'^j is performed in pQCD; in addition 
to the leading order terms, the dominant twist-four term is also included. The 
twist-four term represents multiple scattering in the nuclear medium, and the 
contribution due to this term is a free parameter that can be fixed by compar- 
ison to data. This model was able to describe the HERMES krypton data by 
fitting the HERMES nitrogen data. The modification of the fragmentation 
functions represented by R\j{z,'u) emerges naturally from this description. 
In this approach, a novel quadratic dependence of the energy loss on path 
length is predicted^Hl! originating in the non-Abelian analog of the LPM ef- 
fect in QED 14j. However, the assumption that hadronization occurs outside 
the nucleus will break down at low energies or for more massive hadrons. 

Another recent model invokes deconfinement of nucleons in the nuclear 



5 



medium as a mechanism for the modification of fragmentation functions in 
the nuclear medium [T3|. The basic method follows the rescaling models de- 
veloped in the 1980's. This approach combines partial deconfinement with 
absorption of the formed hadron in the nuclear medium. The nucleonic de- 
confinement effectively extends the range in for which gluon radiation 
takes place, modifying the fragmentation functions accordingly. The model 
successfully describes the EMC data and the HERMES data for pions and 
in nitrogen and krypton, although the prediction for K~ systematically 
deviates from the data. 

A fourth recent model adapts and improves stationary string model (SSM) 
techniques developed in the 1980 's and 1990 's. This model, applied to valence 
quarks, successfully fits the HERMES pion data for nitrogen and krypton, 
and predictions are given for kaons and heavier nuclei ^Hl- As with the pre- 
vious model, the prediction for and K~ attenuations are very similar, 
while the HERMES data show a systematic difference. 

Finally, an approach that emphasizes accurate modeling of the final state 
interactions (FSI) is able to explain the HERMES data without invoking 
modification of fragmentation functions |17||. The struck quark propagates, 
evolves into a hadron, and its interaction with the nuclear medium is given by 
a semi-classical coupled channel transport calculation that has been tested 
in other reactions. These calculations, although in preliminary form, seem 
to have enough flexibility to reproduce the HERMES nitrogen and krypton 
data for charged pions and kaons, and protons and anti-protons. There is 
good qualitative agreement throughout, and quantitative agreement for some 
kinematics. 

3.2 Connections between Models and Observables 

The models highlighted above are generally able to describe the published 
data, in spite of the fact that they rely on quite different physical pictures. 
In order to make good progress in understanding hadron formation in nu- 
clear DIS, it is necessary to discriminate among these models using more 
data that span a wider range of kinematics or that introduce new observ- 
ables. Similarly, more complete calculations are needed. For instance, model 
predictions for the flavor and mass dependence of Rm{z, v) can be tested by 
studying pions, kaons, and protons. In fact, a wide range of particles can be 
considered, as seen in Table ^ In this table, hadrons with cr much greater 
than nuclear dimensions are listed. To the extent that hadronization is the 
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production mechanism, inter-comparison of R\j{z, v) for these particles will 
shed light on such topics as whether semi-classical calculations of FSI suffice, 
or whether a quantum mechanical description is required. Likewise, study 
of moments of transverse momentum may distinguish between models with 
different assumptions about the number of multiple scatters or the shape of 
the emitted gluon momentum distribution. Measurements of R\i at low z 
may be useful, particularly at low energies, to distinguish between pion flux 
entering other reaction channels, as in FSI, from a change in the intrinsic 
fragmentation function, i.e., due to gluon emission. A clean evaluation of 
the dependence of R% will help distinguish between different reaction 
mechanisms. Studies of the angular correlations between two high-energy 
pions may also clarify whether the reaction mechanism is dominated by FSI 
in the nuclear medium or by the fragmentation process. Spin observables, 
such as single-spin asymmetries, may offer indications of whether partonic 
multiple scattering is an important component of the process. Finally, the 
twist-4 pQCD model makes the assumption that hadronization takes place 
outside the nucleus altogether. The point at which this assumption is no 
longer valid is an important unknown that can be addressed by the CLAS 
low-energy data. 

While these are rich and exciting prospects for learning much new physics, 
there are a few experimental and theoretical complications that require fur- 
ther study and clarification. For instance, hadron formation lengths can be 
derived from the data, however, the definition of the formation length de- 
pends on which model is considered. As an example, hadron formation is 
conceived as a multi-step process with multiple time constants in some mod- 
els, but a single-step process in others. Reference frame considerations also 
enter the interpretation of the data; a lower limit on x is required to pre- 
serve the simple picture that a single quark propagates with initial energy 
t/[TH]. Resonances in the residual system may need to be suppressed, which 
reduces the kinematic coverage available. There may be an issue of isolat- 
ing current fragmentation from target fragmentation, although the nuclear 
medium may perhaps suppress the latter. Finally, the degree of validity of 
QCD factorization needs to be evaluated. The last few issues are common to 
all semi-inclusive deep inelastic scattering, and as such are being investigated 
by several physics communities; further progress is therefore to be expected. 
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Table 1 : Final-state hadrons potentially accessible for formation length and trans- 
verse momentum broadening studies in CLAS. The rate estimates were obtained 
from the LEPTO event generator for an 11 GeV incident electron beam. (The cri- 
teria for selection of these particles was that cr should be significantly larger than 
nuclear dimensions, and their decay channels should be measurable by CLAS^"*".) 



hadron 


cr 


mass 


flavor 


detection 


production rate 






(GeV) 


content 


channel 


per Ik DIS events 


7r° 


25 nm 


0.13 


uudd 


77 


1100 


7r+ 


7.8 m 


0.14 


ud 


direct 


1000 


7r~ 


7.8 m 


0.14 


du 


direct 


1000 


V 


0.17 nm 


0.55 


uuddss 


77 


120 




23 fm 


0.78 


uuddss 


7r"'"7r~7rO 


170 


r]' 


0.98 pm 


0.96 


uuddss 




27 




44 fm 


1.0 


uuddss 


K+K- 


0.8 


K+ 


3.7 m 


0.49 


us 


direct 


75 


K- 


3.7 m 


0.49 


us 


direct 


25 




27 mm 


0.50 


ds 




42 


P 


stable 


0.94 


ud 


direct 


1100 


P 


stable 


0.94 


ud 


direct 


3 


A 


79 nini 


1.1 


lids 




72 


A(1520) 


13 fm 


1.5 


uds 


pn^ 




S+ 


24 mm 


1.2 


us 


pn^ 


6 


SO 


22 pm 


1.2 


uds 


A7 


11 




87 mm 


1.3 


us 


Att" 


0.6 




49 mm 


1.3 


ds 


Att" 


0.9 



4 Future Prospects 

Over the next few years it is expected that new datasets will become available 
which will provide further stringent tests of these models. Some HERMES 
data were taken with a 12 GeV beam on both nitrogen and krypton targets; 
27 GcV data on helium and neon are also under analysis. A new program 
at Jefferson Lab will provide even lower energy data at high luminosity on 
targets ranging from carbon to lead. This 5 GeV data, from the CLAS EG2 
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running period, will have sufficient statistics to map out the full multivariate 
character of for several hadron species, and to directly measure transverse 
momentum broadening. 

Other future prospects include the possibility of additional targets in 
HERMES such as xenon, to complete the study of nuclear mass dependence 
at high energy transfer, and a large program which would be possible at 
Jefferson Lab following the energy upgrade. The ultimate moderate-energy 
program will be feasible with CLAS^^, due to the planned upgrade in lu- 
minosity (to 10^^cm~^s~^) and improved particle identification capabilities. 
These improvements will permit mapping out multivariate R'l[ for all the 
hadrons in Table 1 as well as pt broadening for many of these, spanning a 
kinematic range of v from 2 — 9 GeV and 2 — 9 GeV^ in Q"^. 

5 Conclusions 

The study of hadronization from semi-inclusive nuclear DIS is a new topic 
that is drawing increasing interest from several different physics communities. 
The recent wave of interest stimulated by the HERMES measurements has 
resulted in a number of theoretical efforts to describe the data. These efforts, 
based on a variety of physical pictures, demonstrate the exciting possibility 
of gaining new insight into hadronization as well as other aspects of quark 
propagation through nuclei such as quark energy loss. However, more data 
is needed to definitively discriminate among these model approaches. In 
addition, much more theoretical effort is needed to address all of the recent 
data, to reduce the number of assumptions in the models, and to make greater 
contact with the fundamental theory, QCD. 
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